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Abstract 

We propose a novel method to study flavor-changing neutral currents in the e'’'e“ —s- D*^ and 
e'’“e“ —7- B* transitions, tuning the energy of collisions to the mass of the narrow vector 

resonance D*^ or B*. We present a thorough study of both short-distance and long-distance 
contributions to e'’'e“ —)■ D*^ in the Standard Model and investigate possible contributions of new 
physics in the charm sector. This process, albeit very rare, has clear advantages with respect to 
the D^ —)• e'’'e“ decay: the helicity suppression is absent, and a richer set of effective operators can 
be probed. Implications of the same proposal for 5* are also discussed. 


1 



I. INTRODUCTION 


Experimental studies of the flavor-changing neutral currents (FCNC) are among the most 
promising ways to reveal virtual effects of possible new physics (NP) in heavy meson de¬ 
cays. The FCNC transitions have been thoroughly studied in the 6-quark sector, where the 
Standard Model (SM) is seen to dominate the decay amplitudes [HE]- Nonetheless, recent 
hints at anomalies in the exclusive B —decays [3l 0] call for additional studies 
in the 6-flavor FCNCs, preferably, with observables having less hadronic uncertainties than 
in the rare semileptonic decays. The leptonic decays Bs d —)■ £^i~ remain the “cleanest” 
probes, however, only for the axial-vector/pseudoscalar effective operators. Moreover, due 
to helicity suppression, only the £ = /i, r modes are accessible in leptonic decays, leaving the 
detection of the electron modes extremely challenging. 

Still, even if the 6 —)■ transitions are conform with the SM, the situation might 

be different in the charm-quark FCNCs, where ample room for NP effects in the c —>■ u£~^£~ 
transitions is available (see, e.g., 0). Rare decays of the type —)■ £+£ for £ = fi,e have a 
potential to probe a variety of NP scenarios. The studies of these decays are seemingly ap¬ 
pealing because for the local operators mediating these transitions, both in SM and beyond, 
all nonperturbative effects are accumulated in a single parameter, the decay constant. 
However, in SM the long-distance effects dominate the rare decays of charmed mesons. It 
is extremely difficult to estimate these contributions model-independently. Moreover, since 
the initial state in the —)■ e+e“ decay is a pseudoscalar meson, the helicity suppression 

again makes the observation of this process very difficult. This feature persists in many NP 
models as well. In addition, the rare radiative decay —>■ ye+e" has a branching ratio 

that is enhanced by a factor of 0{am‘j^/m1) compared to —>• e+e“, since the additional 
photon lifts the helicity suppression |6]. If the emitted photon is soft, this transition can be 
easily misidentihed as —>■ which further complicates its experimental observation. 

An interesting alternative to — )■ e’''e“ process that is not helicity-suppressed is a 

related decay Zi)*(2007)° —)■ While also probing the FCNC cu —)■ £'^£~ transition, this 

decay is sensitive to the contributions of operators that —)■ £'^£~ cannot be sensitive to. 
Unfortunately, a direct study of the D* —)■ decay is practically impossible, since the 

D* decays strongly or electromagnetically. 

Nevertheless, as we shall argue in this paper, it might be possible to probe the D* — )■ 
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FIG. 1: Probing the cu —>■ e+e vertex with the D*(2007)‘^ resonance production in e+e collisions 


transition experimentally. Assnming time-reversal invariance, it wonld be eqnivalent to 
measnre the corresponding production process e''“e“ —)■ D*, as shown in Fig. [T] In order 
to do so, we propose to run an collider, such as BEPCII [7j, at the center-of-mass 

energy corresponding to the mass of the D* meson. Note that BEPCII already scanned this 
region of energies, achieving the luminosity of about 5 x 10^^ cm“^ s“^ around ^/s = 2 GeV 
[H]. If produced, the D*^ resonance will decay via strong [D*^ —)■ D^tt^) or electromagnetic 
(D*^ _i. interactions with branching fractions of (61.9 ± 2.9)% and (38.1 ± 2.9)% 

respectively.^ 

In the setup discussed in this paper, the D* production process is very rare. However, the 
identihcation of even a single charmed-meson hnal state from the e’''e“ —)■ D*^ —decay 
would provide an unambiguous tag for this flavor-changing production process. Naturally, 
one also needs an adequate quality of the —)■ 27 identihcation and pion-kaon separation 
in the decays in order to reject background processes. 

Our proposal may also be realized in the 6 -quark sector by scanning the region of the 
B*^ resonances at an collider. This will probe the processes B*^ originating 

from the b —)■ s{d)£^i~ quark currents. In fact, studying the transitions involving electrons 
could also shed some light on recent hints at lepton non-universality in 6 — se^e~ versus 
b —)■ |1]. 

Note also that tuning an e’''e“ accelerator to the masses of resonances is not the only 
possibility to access their production. Some sensitivity to these processes could be also 

^ Note that the charged mode D*^ —?> D^tt~ is forbidden by the lack of the available phase space. 
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achieved by studying radiative return events at currently running e’''e“ machines operating 
at their nominal energies. 

The rest of this paper is devoted to a more detailed discussion of this proposal and to 
the relevant theoretical estimates. 

II. e+e- ^ D* RESONANT PRODUCTION. 

Let us consider a generic scattering amplitude of —)■ Dti, and assess the contribution 

of the narrow resonance D* to this process depicted in Fig. Writing this amplitude as a 
matrix element of a generic lepton-quark interaction 

^ic-ff,u){eYe), ( 1 ) 

where only the vector currents are kept for simplicity and an effective scale M and dimen¬ 
sionless coupling A' are introduced, we obtain 

Wl(e+e“-)■ T>°7r°) = (L>°(pD)vr°(p^)|'H|e+(p+)e“(p_)) 

= ^(^°(PD)vr°(p^)|c7^M|0)(0|e7^e|e+(p+)e“(p_)) (2) 

= u(p+)7^u(p_) , 

where the lepton current is factorized out and the hadronic matrix element is expressed via 
the —)■ 7r° vector form factor at s = {po + VttY > (^d + Note that y/s is the 

center-of-mass energy of the e’''e“-collision. Up to an isospin factor, the same vector form 
factor appears in the semileptonic —)• 7r~i^i>£ decay, where s < {mB — m^Y. 

At hrst place, it is the effective coupling in Eq. Q that determines the value of the cross 
section calculated from (|^. Yet, the presence of a narrow resonance in the form factor is 
also crucial. To see that, we isolate the two lowest resonance contributions to the form 
factor, that is, D* and D*' = iA(2600), where the latter, with the mass = 2612 ± 6 
MeV and total width T^./ = 93 ± 14 MeV [^, is the most suitable candidate for the hrst 
radial excitation of iA*-meson. In the resulting decomposition 

n+ ( ') — /£>*°5'D*°D07r0^D*0 _ fD*0'9D*0'D0n0TTT'D*0' _i [ f + Ml 

2(m^*o — s — im^*oro) — s — imB*oiT b*oi) ^ Jbgr ’ 

we expressed the residues of both poles via decay constants of vector resonances and their 
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strong couplings to Dtt, defined, respectively, as 


(4) 

(5) 


{0\u^>^c\D*{p)) = fD*mD*e^{p) , 

{D^{Pd)t^^{p-k)\D*^{p)) = - 5 'D* 0 D 07 r 0 (e(p) ■ p,r) , 

where e^{p) is the D* polarization vector, p = Pd + Pn, and we isolated the background 
contribution [/^o^o(s)] with respect to the two resonances. 

While experimentally only an upper bound for the total width of D^*, Fq < 2.1 MeV, is 
available [U], we can compute the actual value of this width from the measured total width 
of the charged D~^* meson, r_|_ = 83.4 ±1.8 keV [S]. Using the isospin symmetry to relate the 
strong D*Dti couplings and taking into account the phase space correction (see, e.g..[Tn]) 
we obtain: 


Fo 




^ U)°7r°) ± T{D*^ 

r+BD*+^DO^+ f A(m^«o,m^o,m^o) \ ^ Bp^o^po^ \ 

2 + \ Bp*o^po^o J 


~ 60 keV, 


( 6 ) 


where Bp*^f = r(Zl* —>■ f)/TQ is the branching fraction of D* —)■ /, X{x,y,z) is the kine¬ 
matic Kallen function and we employed the data on branching fractions from [9] , neglecting 
small experimental errors. 

The Breit-Wigner form used in Eq. ([^ is certainly applicable for such a narrow resonance 
as D*. We adopt the same approximation for the broad resonance D*\ which is sufficient 
for an order-of-magnitude estimate. At s = m|,* the magnitude of the excited state con¬ 
tribution to the form factor is suppressed with respect to the Zl*°-pole term by the factor, 
approximately. 


fpo*/ yp*'opo^o'fnp*oi 


iTo 

foo* 90*0 DO 170^0*0 


2 A — iT 


5.0 ■ 10 


-5 


(7) 


where A = m^,o' — mD* ~ 600 MeV and the 0{A/mp*) terms are neglected. In the above, 
we assume that the strong couplings and decay constants of the radially excited and ground 
D* states are in the same ballpark, so that the first factor in this ratio is of 0(1). In fact, 
the ratio of decay constants in Eq. (|^ is less than one, as the QCD sum rule predictions 
indicate m . Our estimate also implies that if one produces a resonance with a width of 
0(100 — 150) MeV, typical for the light vector mesons such as p or K*, there is no relative 
gain in the resonance cross section. We emphasize that the very small width of D*, driving 
up the cross section, essentially originates from the miniscule phase space of its single strong 
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decay mode. In this situation, the radiative mode of D* becomes equally important and the 
FCNC transition acquires a tiny but non-negligible branching fraction. 

Having in mind a strong suppression of all other than D* contributions to Eq. (|^ at 
i/s ~ m^)*, the cross section of e’'“e“ —?■ Dtt can be written in a standard resonance form: 


a 


(e+e ^ = aD*(s) = 


127r 


B 


m 


D" 


e~ B£)* 


mlSl 


D* 


(s - 


( 8 ) 


To exploit the D*-resonance enhancement around ^/s = rriD*, described by Eq. ([^, we 
tacitly assume that an appropriate tuning of the electron-positron accelerator beams can be 
performed, so that their energy resolution is smaller than the spread of the resonance cross 
section. Then we can simply use Eq. ([^ at s = m'jj,, yielding 

127r 


aD*{mjy,) = B 


D*—^e”*"e Bd*^d-j^ 2 


m 


~ HB*^e+e-(2.26 X 10^) nb. 


(9) 


D* 


Let us recall that the total cross section a{e~^e —)■ hadrons) is about 50 nb at y/s = 2.0 

GeV [g. 

The expected number of e’''e“ —)■ Dti events at i/s = tud* is given by the product 


Nd* = e / Ldt 


( 10 ) 


where e and J Ldt are the detection efficiency and time-integrated luminosity, respectively. 
The condition > 1 leads to a lower bound on the D* —)■ branching fraction that 
still allows one to detect the process —?■ —)■ Htt: 


B 


D*->e+e“ 


> 



m 


X 


D* 


127r B 




( 11 ) 


For example, an average luminosity at the level oi L k, 1.0 x 10^^ with a 

’’Snowmass year” (~ lO^s) of running at the D* resonance yields J Ldt = 1.0 fb~^. 

Under these conditions, the single-event sensitivity implied by the bound ( [IT| ) means that 
branching fractions 

BD*^e+e- > 4 X 10"^^ (12) 


could in principle be probed. While the above bound is a very crude estimate of the possible 
sensitivity, as we did not take into account detection efficiency, assuming e = 1, it can serve 
as a useful criterion. In order to see if and when this bound can be approached, let us hrst 
calculate the branching fraction Bu^^e+e- in fhe SM. We need this quantity as it is the key 
parameter determining the —)■ D* cross section. 
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III. D* e+e- TRANSITION 


The most general expression for the D* —)■ e'''e decay amplitude can be written as 


A{D*^e^e )=u{p-,S-) 


C 

+-(p+ - p_) 


TTl £)* 


D 


+- {p+-P-)Alb 


ruD 


v{p+,s+) e^(p), 


( 13 ) 


where A, B, C, and D are dimensionless constants which absorb the underlying effective 
quark-lepton interaction and the vacuum —)■ D* hadronic matrix elements (e.g., the D* decay 
constant). Note that the simplified effective interaction ([^ used in the previous section 
corresponds to A = {X'/M‘^)f£,*m£)* and all other constants put to zero. The amplitude 


(13) leads to the branching fraction 

rriD* 


B 


D*^e+e- — 


127rrr 


{\Af+\Bn+^{\cf+\Dn 


(14) 


where we neglected the mass of the electron. 

It is straightforward to compute the SD part of the D* —>■ amplitude triggered by a 

generic effective Hamiltonian, not necessary containing only the SM operators. The decay 
amplitude is 

(e+e-ineslD*) = GY,c^{fi){e+e-\Q,\D*)\,, (15) 

i 

where G is a constant with the dimension of inverse squared mass that sets the scale for the 
operators. For example, G = AGp/\/2 in the SM or G = 1/A^ for new physics. In the above, 
Qi are the effective operators of dimension six, and q are the corresponding (dimensionless) 
Wilson coefficients. The most general set of ten local operators producing the cu —)■ £'^i~ 
transitions reads 


Ql = {Gl^,£L) Qa = {£r£l) {urCl), 

Q2 = (uR'y^CR), Qs = (Ir<j^i,£l) {uRa^’^CR), 

Q3 = (£l£r) {urCl), 


(16) 


with hve additional operators Qq, ... ,Qio obtained respectively from those in Eq. (16) by 


the substitutions L ^ R and R ^ L. We will be using the set of operators Qj=i,.,,,io to 
assess possible contributions of any generic new physics model to D* —)■ 
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The operators with (pseudo)scalar quark currents Q 3 - 4 , and Qs-g do not contribute to 
this process. On the other hand, contrary to the —)■ decay, not only the vector 

operators Q 1-2 and Qq-j but also the tensor operators Q 5 and Qiq could, in principle, 
contribute to D* —)■ e+e“. However, in most of NP models, these effective operators are 
absent, motivating us to neglect their contributions. This is equivalent to setting C = D = 0 


in Eq. (13). In terms of the Wilson coefficients q the constants A and B are 

G 

A = —(ci + C 2 + Cg + Cj ), 


G 


B = - — (Ci + C 2 


C 6 


C 7 ). (17) 

The effective operators in SM (see [T3] for the definition) mediating the D* —)■ e+e“ transition 


in SM are easily matched to the operator set in Eq. (16), 


Og = 


IGtt^ 


+ Qr 


Oio — 


IGtt^ 


Q] 


(18) 


The expressions for the corresponding Wilson coefficients Gg and G^g, where a superscript c 
indicates that CKM matrix elements are included in their definitions, can also be found in 

ra- 

In addition to O 9 and Oio, in SM the D* —)■ amplitude receives a contribution from 

the magnetic dipole operator O 7 coupled to the leptons via virtual photon. The correspond¬ 
ing part of the effective Hamiltonian reads 


rr(77) _ „c,eff 


rTTlc UlO-^’^CrF, 


flU 


(19) 


^/2 VlGvr^ 

A calculation of this contribution requires the knowledge of the tensor (transverse) decay 
constant of the D*\ 

(0|mct^^c|T)*(p)) = ifl. (eV . (20) 

In the absence of the estimate of /^*, we rely on the properties of light vector mesons for 
which the vector (longitudinal) and tensor decay constants are in the same ballpark and 
simply assume that = f^*. One can now calculate the SD contribution to D* —>■ 
in the SM in terms of effective constants: 

me /J 




O 9 ’"® + 2 


D* QC,eS 


^D* fo 


( 21 ) 











which results in the SD contribution to the branching fraction: 


B 


a 


2r<2 




ObTr^Fo 




D* 


fT 

QC,eS 2 J D* 


ttld* Id* 


+ |C^: 


101 


( 22 ) 


Adopting rric = 1.3 GeV, we use for the Wilson coefficient^, Cg(fi = rric) = 0.198|14\V);b|, 
neglect CIq{^ = rric) and employ the results of the two-loop calculation [H] for the remaining 
effective coefficient, = rric) = —0.0025. We also use the central value of the QCD 

sum rule estimate fa* ~ 242 MeV [15]. Substituting all input paramaters in Eq. (22), we 
find 


B 


SD 

D*^e+e~ 


2.0 X 10“^®. 


(23) 


As expected, this number is extremely small, several orders of magnitude below the lowest 


accessible branching fraction Eq. (12) and thus beyond any realistic experimental setup. 
Still, it is instructive to remind the reader that the short-distance width of the similar decay 
of the pseudoscalar is many orders of magnitude smaller: ~ 10“^^, whereas 

~ 10-^® (see, e.g., HH). 


IV. LONG-DISTANCE CONTRIBUTIONS 


Generally, in rare charm decays, a significant enhancement of the decay rate is expected 
in SM due to LD contributions, generated by the four-quark weak interaction combined with 
the emission of the e+e“-pair via virtual photon. It is very difficult to reliably estimate these 
contributions in U —)■ decay because the two-photon intermediate state overlaps with 

long-distance hadronic interactions. 

To investigate the case of D* —)■ decay, we isolate the relevant AG = 1, single 

Gabibbo-suppressed transitions in the effective Hamiltonian T-Lu, of the SM, representing it 
in a form of the two four-quark operators: 


= 


4G; 


E 


q=d,s 




c: 


'c{q )' 


Nc 


{<lLluqL){uLlyCL) 


+2cy>(5i7..r“<;L)(t<L7.r“ci) 


( 24 ) 


^ We thank Dirk Seidel for providing us this coefficient calculated at NLL order. 
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where T“ are the color-octet matrices. Note that we again included relevant CKM matrix 
elements into the dehnition of and 

The LD contribution to the D* —>■ e’''e“ decay amplitude is given by the matrix element 


(e+e \n^,\D*{p)) = -e m(p_, s_) 7 ^n(p+, s+) 


p2 


(25) 


9 2 


where we factorized out the lepton current, the photon propagator and dehned the hadronic 
matrix element in a form of the correlation function 


E,(p") = * / d‘xe‘’‘‘{0\T {jy(i)W„(0)} \D’{p )), 


(26) 


where j®™' = 'Yl,q=udscQ<i^l^d^ electromagnetic (e.m.) quark current. The quark-level 

diagrams corresponding to this amplitude are depicted in Fig. where we distinguish two 
main topologies corresponding to the virtual photon interacting with (a) the q = d,s quarks 
(upper panel) and (b) with the u quark (lower panel). Both contributions develop imaginary 
parts due to the intermediate light-hadron states. We neglect the e.m. interaction with the 
heavy c-quark. 

To estimate the LD contribution corresponding to the diagram Fig. we adopt factor¬ 
ization approximation, that is, neglect the gluon exchanges between qq and he helds in 


Hence, we retain in Eq.(26) only the product of color-neutral vector currents and obtain: 


sjf'(p=) = ^ E <21 ^ /<iVe-'-(0|r{57,9(1) 97,5(0)} |0>| 

^ q=d,s \ ^ 

X {0\uYc\D*ip)), (27) 

where we recognize the quantity in the curly braces as the polarization tensor for the 
g-flavored vector current. Substituting its decomposition 


njfj(p) = {-gfduv^ + v^dVv) (p^) (28) 

in Eq. (27) and parametrizing the matrix element of the uc current according to Eq. (|^, we 


hnally obtain the LD transition amplitude in the form (13) with 

= 4iraQ,^/D.mD. (^C{ + (n^)?'^) - n“>(p^)) 


(29) 


9 9 
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FIG. 2: Long-distance contributions to e'^e~ —?• D*^ caused by the virtual photon interaction (a) 
with q = d,s quark in factorizable approximation, forming the polarization operator; (b) with u 
quark. 


and = 0. In the above expression, we slightly modihed the Wilson coefficients, so 

that C ^2 = taking into acconnt that the CKM factors for the s- and d-qnark 

parts of T-Lw are approximately eqnal and have opposite signs. As expected, in the flavor 
SU (3) limit the whole LD contribntion vanishes, reflecting the GIM cancellation. 

Note that the factorizable qnark-gluon effects are implicitly retained in the two separate 
hadronic qnantities, the decay constant and the polarization operator. Nonfactorizable QCD 
corrections will involve contribntions of the color-octet and axial-vector qnark cnrrents. We 
expect that these corrections are snppressed either by asijric) or by the powers of Kqcd/'^o 
since the characteristic momenta flowing throngh the qnark loop in the diagram of Figj^ 
are of order of ~ ^c- For the same reason, one may also argne that the mixing between 
ss and dd loops in the polarization operator is snppressed, which is also in line with the OZI 
snppression valid in the vector-meson channel. To snbstantiate the expected snppression of 
higher-order effects, one has to calcnlate the nonfactorizable mnlti-loop diagrams explicitly, 
a technically difficnlt task we postpone to fnture studies. 


Returning to the LD contribution (29) we estimate hrst the polarization operators. Note 


11 








that satisfies a (once-subtracted) dispersion relation 




p 


ds- 




12ti‘^Q'1 Jq s{s - p^ - ie) ' 


(30) 


where R^^\s) is the normalized e+e cross section to hadrons initiated by the quark fiavor 
q = s,d, so that below the charm-anticharm threshold 


R{s) 


cr(e+e —)■ hadrons) 
cr(e+e“ —)■ 




(31) 


In principle, Rq{s) could be extracted from the experimental data. In practice, however, it 
is hard to disentangle hadronic states generated by the d and s-quark currents. Therefore, 
for the sake of estimate, we employ the parametrization of Rq{s) stemming from the QCD 
sum rule analysis UTl of the M, d- and s-quark correlation functions and based on the quark- 
hadron duality: 


R^^\s) 






E 

V=p^,U) 




a. 




TT 



(32) 

(33) 


where the decay constants of the light vector mesons are defined via matrix elements of the 
electromagnetic current: 


{0\jr\V) = nvmvfve'(yy (34) 

Here the coefficients Kp = l/\/2, = 1/(3 a/2) and = —1/3 reflect the valence quark 

content of these mesons. Furthermore, we take into account the total widths of the vector 


mesons, replacing delta-functions in Eqs. (32 33) with the Breit-Wigner approximations. 


To calculate the difference of the polarization operators entering Eq. (29), we use fp = 220 
MeV, ft^ = 197 MeV, and = 228 MeV, obtained from the experimental values of the p, u 
and 0 leptonic widths [H], We also choose the effective thresholds [H]: Sq = 1.5 GeV^ 
and Sq = 1.95 GeV^. To finalize our estimate of the LD amplitude, we take the relevant 
combinations of GKM parameters pckm = Vus^cs — —V*^Vcd = A(1 — A^/2), where A = 
0.22537 [9] is the Wolfenstein parameter, and use the Wilson coefficients ai p = me = 1.3 
GeV at LO(NLO): Cl/pcKM = —0.53(—0.41) and C^/pckm = 1.28(1.21). 

If we assume that the LD contribution of the annihilation type calculated above saturates 
alone the decay amplitude, the resulting branching fraction of the D* —)■ decay, turns 
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(36) 


out quite sensitive to the mutual cancellation of the Wilson coefficients C\ and C 2 / 3 : 

„LD.A ^ / 4.7 X 10-» (NLO) 

\ 5.7 X 10-‘* (LO) 

The cancellation is numerically less pronounced in LO, in which case the above branching 
fraction grows by almost two orders of magnitude with respect to the NLO case, becoming 


substantially larger than the SD width (23) 


One could also try to estimate the difference of polarization operators in Eq. (29) by 


subtracting the d- and s-quark loop diagrams from each other. This is equivalent to replacing 


by the parton spectral density (the second term in brackets in Eqs. (32 33)) and 
integrating from to inhnity. In this case, to account for a correct normalization of 
the effective operators in and to avoid infrared unstable terms of 0 {ln{md/ms)) in the 
difference of the polarization operators, one has to add a (scheme dependent) constant term 
to the dispersion integral, which contains In The resulting difference of the loop 

functions is then proportional to (m^ — m^)//i^ as it should be in the GIM cancellation.^ 


The simple loop estimate yields, instead of Eq. ( |35| ), 


- rs> 


2.0 X 10-20 (NLO) 
2.4 X 10-1® (LO) 


(36) 


in the same ballpark as the estimate (35). 


Turning to the LD contribution represented by the diagram in Fig. |^b, we present a 
rongh estimate of its amplitnde, taking the imaginary part and accounting only for the 
lowest TT+TT- and K^K~ intermediate states. Other contribntions, including mnltiparticle 
ones, are possible. Those can in principle give even larger amplitudes due to phase space 
suppression of the S'17(3)-related intermediate states [22] . 

According to the unitarity condition, the imaginary part of the D* —)■ e+e- amplitude 
generated by these particular intermediate states is expressed via the e’''e- —)■ P^P~ e.m. 
amplitnde which contains the {0\QuU''yf^u\P~^P~) form factor taken at the timelike momen¬ 
tum transfer The e.m. amplitude is multiplied by the D* —)■ P^P~ nonleptonic 

amplitude (P = vr, K) and integrated over the P^P~ phase space. 


^ In the hadronic representation of the polarization operator used above, the same cancellation can be traced 
in the difference of QCD sum rules for and /| and stems from the difference of loop diagrams in 
the perturbative part and, in addition, from the 0{ms{ss) — md{dd)) terms originating from the vacuum 
condensate contributions. 
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Applying the naive factorization approximation for the nonleptonic amplitude, we intro¬ 
duce the relevant hadronic matrix element 


Pt(P (P-)\lT'1sc\D'(p)) = 2i{e{p) ■ , 


(37) 


which depends on the one particular form factor of the D* ^ P transition, (we dehne these 
form factors analogous to the well familiar D ^ V form factors). Using also the decay 
constant of the light pseudoscalar meson: 


(F+(p+)|M7^75g|0) = -ip%fp, 
we obtain the LD contribution in the following form: 


(38) 


X 


1 - 


fKf^K 


D* 


(39) 


where jSp = For an accurate numerical estimate we need to calculate the 

form factors of D* ^ P transition, which is beyond our scope, demanding, e.g., a dedicated 
application of QCD light-cone sum rules. For an order-of-magnitude estimate we assume that 
these form factors are in the same ballpark as the (correspondingly normalized) D ^ 7 t,K 
form factors and take A^*’^(m^) ~ 1. The pion e.m. form factor F®™'(m|)*) ~ 0.28 is known 
from a measurement [18] in the time-like momentum region. For the SU(3)-violating ratio 
of the e.m. and heavy-light form factors we conservatively assume that they can vary within 
±30% 

Af“('mj()Fp'{ml.) ^ , , „ , ,. 

A^’’(ml)Fr(ml.) ' ' ' ' 

After that we obtain the lower limit on the branching fraction assuming that the LD con¬ 


tribution (39) is dominant: 


f^{LD,b) 

‘^D*^e+e 


_ > (0.1 -5.0) X 10 
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(41) 


where the broad interval mainly reflects the variation within the limits assumed in 
whereas switching from NLO to LO Wilson coefficients is a minor effect in this case. 

A more complete analysis of LD effects including possible interference between contri¬ 
butions of the two types considered above as well as a more complete set of intermediate 


14 







states is an interesting task we postpone for the future. The estimates presented above 
demonstrate possible approaches which could be developed further. 


The branching fractions presented in Eqs.(35), (36) and (41), can achieve the level of 
10“^®, that is, not very much larger than the typical values (23) in the presence of only 


SD effects. Altogether, the probability of e^e —)■ D*^ in SM remains considerably smaller 


than the estimated minimal level (12) reachable (optimistically) by experiment. Therefore, if 


observed or at least constrained at the level of Eq. (12), the e^e —)■ D* events will dehnitely 


have no SM background. While this is basically true for a majority of FCNC transitions in 
the charm sector, such as H —)■ or semileptonic decays, the —)■ D* transition has 
a clear advantage of having a relatively moderate LD background. 


V. D* 


e^e 


TRANSITIONS AND NEW PHYSICS. 


It is interesting to estimate what possible NP scale could be probed by the processes 


discussed in this paper. As follows from Eq. (16), there are ten possible operators that 


parameterize any NP contribution to any cu — )■ e’''e“ process. Some of those operators can 
be probed both in —)■ e’''e“ and D^* —)■ e^e“ transitions. Some can only be reached 

in D^* —)■ or (more challenging from hadronic point of view) rare semileptonic charm 

decays D —)■ Me^e~. It would be interesting to compare the possible reach of and D^* 
decays. 

Assuming that a NP contribution is dominated by a single operator, it is easy to see from 


Eqs. (14) and (17) that 


A 


1 C‘^ 


37r 32Pn B 


D* 


R + P. 


1/4 


(42) 


where C = |cj| for i = 1, 2, 6, or 7. Employing the upper bound (12) we hnd that observation 


of a single event in a “Snowmass year” of running would probe NP scales of the order of 
A ~ 2.7 TeV provided that C ~ 1. Let us compare this to the NP reach of —)■ 
decay for the same operators. Using [12], 

1/4 


A 




(43) 


\ d,2'KT£) Bo^e+e-, 

where Pd is the total width of the meson, and the current experimental bound, 
Bo^e+e- = 7.9 X 10“®, we find that only scales A ~ 200 GeV are currently probed by 
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D —)■ e~^e~ decay. It is the presence of the lepton mass factor that severely limits the NP 
scale sensitivity in this process. 

To exemplify this discussion, let us consider two particular models of NP to see how well 
they can be probed in D* —)■ e''“e“ transition. The selected models are by no means unique. 
The chosen examples simply illustrate the differences in sensitivities between —)■ 

and D* —)■ e+e“. 

R-parity violating (RPV) SUSY. R-parity violating SUSY models can be probed in FCNC 
charm decays [ni[i3]. The relevant part of the superpotential can be written as 


= A'„, 




(44) 


where the coupling parameters are defined such that i denotes a generation number for 
leptons or sleptons, j ~ for up-type quarks, and k - for down-type quarks or squarks. As 


can be seen from Eq. (44), the relevant FCNC quark transition c + u —>■ e+e is mediated by 


a tree-level d-squark exchange. Since its mass is much larger than the energy scale at which 
D* decay takes place, it can be integrated out, resulting in the effective Lagrangian [El 13], 

AC. Ah 


rnp _ 


d-R 


Qi 


(45) 


which implies that Gci = ACfcACfc/(2mjj,) in Eq. (15). Extracting the coefficients A and B 


and using Eq. 0 leads to the branching fraction 

1 m%fl. 


tqnp 

^D*^e+e- 


3847r mJ,Fo 

R 




(46) 


Contrary to the case of —)■ , no helicity suppression (factors of mj/m 


D, 


IS seen m 


Eq. (46), which in principle makes this process more sensitive to the NP parameters. This 


feature exists for any NP model that is represented by U ± A interactions. Numerically, 
taking updated bounds on AC^ACfc from [19], conservatively (see also [20]), 


AC.A( 




< 3.83 X 10 


-3 


"^4 


(47) 


we estimate that 


JOO GeV_ 

< 1.7 X 10“^^, implying that some prospects exist for improve¬ 


ment on this bound using the described process. 
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Models with Z'-mediated gauge interactions. Another interesting and representative 
model that we would like to consider here is a model with flavor-changing Z'-mediated 
interactions. In general, 


I^Z' — —g'z’l^L'^u^LZ'^ — g'zi2^R'lti^RZ'^ 

- gTnULl,CLZ'>^ - g7'2ULl,CLZ'^^. 


For ruz' ^ rno the Lagrangian in Eq. (48) leads to 

1 


r^' - 


Ml I 


9z'i9z'iQi + 9z'i9z'2Q‘2 + 9z'29z'2Qs + 9z'29z'lQ7 


(48) 


(49) 


Again, identifying the Wilson coefficients q from Eq. (49) and computing A and B leads to 
the following branching fraction, 

1 




D*^e+e~ 


127r M^/Fq 


4 -p \9T'1 + fi'z'21 ( Ifi'z'll + lfi'z'21 


(50) 


As with our previous example, Eq. (50) does not exhibit helicity suppression of the rate. 


Most importantly, is non-zero for purely vectorial interactions of the Z', which 

will be realized if, for example, g'^/^ = g^z' 2 - This is contrary to —>■ decay rate, where 

such vectorial couplings are forbidden by vector current conservation [12]. There are five 
parameters that describe generic Z' interactions with quarks and leptons, g^z’ 2 i 9'z'n 
g'zi 2 i Mz'. To assess the sensitivity of the e''“e“ —)■ D* production mechanism to Z' 
models numerically, let us make two simplifying assumptions. First, let us assume that Z' 
only couples to left-handed quarks,^ which would mean that g^z '2 — 0- Second, let us assume 
that the Z' has SM-like diagonal couplings to leptons. 


9z'i — 


-- +sin^ 6 


cosOiv \ 2 ' cos^vi/ 

where g is the SM SU{2) gauge coupling. The branching fraction would then only depend 
on the combination g'^z'i/Ml^ 


, _ g sin^ Ow 

9z'2 — 


(51) 


B^' 

Or* 


\/2Gf 3 p 2 fl 


SttFc 




D* 


Mi, Ml \A 


— sin 6*w -I- 2 sin 9 


w 


Taking the constraint Mz' / ■\/g‘^z'i > S.7 x 10^ GeV from —)■ /i+p [T2| yields 

< 2.5 X10-“, 

which is far above the SM predictions for this rate. 


(52) 


(53) 


^ Equivalently, we could have assumed that Z' only couples to the right-handed currents. Then g'^Z'i = 0 
and constraints would be obtained for g‘Z' 2 - 
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VI. IMPLICATIONS FOR Bt FCNC DECAYS. 


Similarly to —)■ e'^e~, the Bs(^d) e+e“ decay is helicity-suppressed. Hence, it might 

be interesting to see if a e'''e“ production process can be used to probe the t e'''e“ 

transitions. Hereafter we will concentrate on the B* resonant production, the corresponding 
process with B^ is very similar, but CKM suppressed and has therefore less chances to be 
detected. There are two important differences between the beauty and the charm case. First 
of all, the strong decays of the B* are kinematically forbidden, and the dominant channel is 
the radiative one, B^* —)■ B^j. This feature is welcome because it leads to a smaller total 
width than for D*. However, additional challenges might emerge for triggering the hnal 
state of the produced B'^^y Second, in the SM the B* —>■ decay is dominated by the 

SD contributions, stemming from the well dehned effective Hamiltonian. Hence one may 
expect a reasonable accuracy in predicting the decay rate, also because the only hadronic 
parameter involved is the B* decay constant. On the other hand, one also has to assess the 
possible LD contributions, e.g., the effect similar to the one shown in Fig. but with the 
c-quark polarization operator. 

For e’''e“ —)■ H*, we use the resonant cross section 


127r 



(54) 


and hnd that with the branching ratios 


B(B; ^ e+e") > 2.0 x 10"'^, 


(55) 


the resonance production can be observed with at least one event. Here we again assume 
one year running at the B* resonance energy, with the luminosity of ~ 10^^ cm“^ s“^ and 
e ~ 1 detection efficiency. While no e’''e“ collider is currently operating at this energy, 
future upgrades of the existing e+e“ machines, including also the radiative return setup, 
might make this energy region available for experimental studies. 

Before calculating the branching fraction of B* —)■ e^e“, it is important to £x the to¬ 
tal width of B* which practically coincides with the width of the single flavor-conserving 
radiative decay: 



(56) 
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Here we use the definition of the B*Bs'y coupling: 


{B,{ph{q)\BHp + q)} = v^iro Sb.-b.t , 


(57) 


The dominant contribution to the H* —)■ H'y {H = D,B) transitions stems from the long¬ 
distance photon emission off the light quark in the heavy meson. The analyses of these 
couplings in terms of heavy-hadron ChPT and related approaches [2S1121] allow for a sim¬ 
plified parametrization of the coupling ® : 


Qq I Qq 

gH*H'r — -^-, 

rriH* Pg 


(58) 


where Qqiq) is the charge factor of the heavy(light) quark Q = c,b {q = u, d, s) in H^*\ and 
fiq is a nonperturbative parameter, which does not scale with the heavy mass. Numerically, 
this relation describes well the two experimentally measured D'y widths, if iJ,u,d — 

420 — 430 MeV is taken. Using the same value of qLu,d for the gs^oBO'y coupling, we obtain 
T(B*^ —)■ B^j) ~ 0.2 keV, in the ballpark of the estimates obtained in [2HI21]- To account for 
the S'f/(3)-fiavor symmetry violation, we adopt the model of [21| where the photon emission 
from the light-quark is described via vector-meson dominance, so that qig — 


We obtain then from Eq. (56) 


Tf, ~ 0.07 keV. 


(59) 


Our conclusion is that the i?*-resonance is considerably narrower than D*. 

To estimate the probability of B* —)■ e’''e“ in SM, we employ the relevant SD part of the 
effective Hamiltonian for the b —>■ si^i~ transitions: 


Tdes — 


4Gp ^ 

■^VaV;, Y. CtOi + h.a, 

V ^ .,-_7 Q in 


( 60 ) 


involving the operators: 


Of — 

Og = 

Oio = 


2=7,9,10 


CUlf, _ 


1677^ 


1677^ 

1677^ 


{lipid) 


(61) 


® The H* —)■ Hy couplings are also estimated [2^ from QCD sum rules, relating them to the electromagnetic 
susceptibility of QCD vacuum, the latter parameter is however known with a rather large nncertainty. 
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Using Eq. (22) with the obvious substitutions of charm meson parameters by the beauty 
ones, we obtain for the branching fraction: 




o?G% 


|r,v, 


is \ 


^ Cg + 2 


fB* ^eff 
t ^7 

ms* /b* 


2 



(62) 


To rehne our estimate, we also should add the nonlocal contribution generated by the 
combination of the current-current weak operator Oi ^2 and the quark e.m. current. One 
of the important effects is the lepton pair emitted from the intermediate cc pair, described 
by the diagram similar to Figj^. In the factorizable approximation, this contribution is 
estimated in a full analogy with the LD effect for the charmed vector-meson leptonic decay, 
i.e., we can use the same expression but replace the difference of d and s polarization 
operators with the single c-quark one: 




AiraQc 


Gf 

V2 


ycbV*gf b^TTIb* 





7h(ln^ + l)+n'yp 





(63) 


where the constant term reflects the correct renormalization of the effective operators. The 
corresponding u-quark loop effect is CKM suppressed and therefore neglected. The c-quark 


polarization operator in Eq. (63) at far above the charm-anticharm threshold 

can be estimated using the dispersion representation for the simple c-quark loop diagram 
(“global duality” approximation) 


OO 

- 6 / 


ds 


4m? 


s[s — — le) 




(64) 


Note that the amplitude (63) can be cast in terms of an effective process-dependent addition 


ACg to the coefficient Cg in the branching fraction (62): 


AU, 


Bt^e+e- 




(65) 


in full analogy with the analysis of the nonlocal charm-loop effect in the semileptonic decay, 
such as B —)■ Kii (see, e.g., ffl). At rric = 1.3 GeV, /i = 4.5 GeV and at the Wilson 
coefficients taken at the same scale: Ci{fi) = —0.255, G 2 (/i) = 1.11, the numerical calculation 
yields: 
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revealing a very small correction to the short-distance coefficient Cg. We therefore skip the 
other LD effects, e.g., the one with the e.m. interaction of the s-quark in 5*. 


Finally, to estimate the branching fraction (62) numerically, we replace Cg —)■ Cg -|- ACg, 


use the numerical values for the relevant Wilson coefficients at the scale /i = 4.5 GeV: 

= —0.316, C'g(/i) = 4.293, Gio(/i) = —4.493, and employ the QCD sum rule estimate 
[T5] for the B* decay constant fs* — 250 MeV. We arrive at the following prediction for the 
branching fraction: 


— 0.98 X 10 


-11 


(67) 


This estimate implies that already within the SM, one could expect several events of the 
type —)■ Bgj to be observed. The signature of the hnal state is a combination 

of monochromatic low-energy photon and a flavor-violating Bg meson. 


VII. CONCLUSION. 

We argued that the rare leptonic decays of heavy mesons H* = D*,B* can be probed in 
the reverse process of the e’''e“ —?■ H* production, provided that the beams of e''“e“ collider 
are tuned in resonance with niH*- 

We calculated relevant transition rates for charm and beauty modes. In the case of 
charmed mesons we paid particular attention to the LD effects, which are calculable and 
found them to exceed the typical SD rate by at most one order of magnitude. We also studied 
several examples of NP scenarios and considered similar production effects for the beauty 
mesons. More efforts can be invested in improving the accuracy of the estimates presented 
in this paper, by calculating the QCD corrections to the LD effects and by extending the 
set of NP scenarios sensitive to the processes we considered here. 

It would be interesting to note that similar single-charm (single-^) hnal states can be 
produced in non-leptonic weak decays of heavy quarkonium states, such as J/'ip —)■ Dvr. 
Since heavy charmonium states lie reasonably far away from the energy region discussed in 
this paper, these transitions will not be producing any backgrounds for e+e“ D* ^ Dtt, 
but can be used to study experimental systematics associated with such hnal state. In 
addition, these transitions are interesting on their own and will be discussed elsewhere. 

Although the experimental setup suggested here might look futuristic, we are convinced 
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that the continuous progress of collider and detector technique will make the tasks suggested 
in this paper real. We hope that our hrst exploratory estimates will simulate dedicated 
experimental studies of the heavy vector-meson production in electron-positron collisions. 

Note added 

While we were hnishing this paper, the work IZ7I appeared where similar considerations for 
i?*-meson were presented. 
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